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1 INTRODUCTION 



Recent studies in observational cosmology, using all the 
available high quality cosmological data (Type la super- 
novae, cosmic microwave background, baryonic acous- 
tic oscillations, etc), converge to an emerging "standard 
model", which is flat and it is described by the Fried- 
mann equation: H^{a) = 87rG[pm(a) + Pq{ci)] /3, with 
a{t) the scale factor of the universe, pm(a) the density 
corresponding to the sum of baryonic and cold dark 
matter and an extra component pQ{a) with negative 
pressure called dark energy and needed to explain the 
observed accelerated cosmic expansion (eg., Davis et al. 
2007; Kowalski et al. 2008; Komatsu et al. 2009; Hicken 
et al. 2009 and references therein). 

The nature of the dark energy is currently one of 
the most fundamental and difficult puzzles in physics 
and cosmology. Indeed, during the last decade there has 
been an intense theoretical debate among cosmologists 
regarding the nature of the exotic "dark energy" . Due 
to the absence of a physically well-motivated fundamen- 
tal theory, various candidates have been proposed in the 
literature, among which a cosmological constant (con- 
stant vacuum), a time varying vacuum quintessence, 
fc— essence, vector fields, phantom, tachyons, Chaplygin 
gas and the list goes on (eg., Ozer & Taha 1987; Wein- 
berg 1989; Wetterich 1994; CaldeweU, Dave & Stein- 
hardt 1998; Brax & Martin 1999; Peebles & Ratra 2003; 
Brookfield et al. 2006; Boehmer & Harko 2007 and refer- 
ences therein). The simplest type of dark energy corre- 
sponds to a scalar field having a self-interaction poten- 
tial V{(f)), with the field energy density decreasing with 
a slower rate than the matter energy density (dubbed 
also "quintessence", eg. Peebles & Ratra 2003 and ref- 
erences therein), and the dark energy component be- 
ing described by an equation of state pQ = wpg with 
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w< —1/3. Note, that a redshift dependence of the equa- 
tion of state parameter is also possible but its present 
functional form is phenomenologically based (see Cheva- 
lier & Polarski 2001; Linder 2003). A particular case of 
"dark energy" is the traditional cosmological constant 
(A) model (corresponding to w= —1), which appears to 
be supported by the combined analysis of the recent rel- 
evant observational data (eg. Komatsu et al. 2009 and 
references therein). 

It has been shown that the application of the cor- 
relation function analysis on samples of high redshift 
galaxies or X-ray selected AGN can be used as a useful 
tool for cosmological studies (eg. Matsubara 2004; Basi- 
lakos & Plionis 2005; 2006). The scope of the present 
study is along the same lines, ie., to place constraints 
on the (r2m,w) parameter space of spatially flat cos- 
mological models using a single cosmologically relevant 
experiment, ie., that of the recently derived clustering 
properties of the XMM-Newton soft (0.5-2 keV) X-ray 
point sources (Ebrero et al. 2009a). 



2 OBSERVED AND PREDICTED 
CORRELATIONS 

2.1 X-ray AGN Correlations 

Recently, Ebrero et al. (2009a) derived the angular cor- 
relation function of the soft (0.5-2 keV) X-ray sources 
using 1063 XMM-Afewton observations at high galac- 
tic latitudes (hereafter 2XMM). A full description of 
the data reduction, source detection and flux estima- 
tion are presented in Mateos et al. (2008). In brief, the 
survey contains ~ 30000 point sources within an effec- 
tive area of ~ 125.5 deg^ (for an effective ffux-limit 
of fx > 1.4 X 10~^^ erg cm~^ s~^ ). Also, Ebrero et 
al. (2009a) presents the details regarding the angular 
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correlation function estimation, the various biases that 
should be taken into account (the amplification bias and 
integral constraint), the survey luminosity and selection 
functions as well as issues related to possible non-AGN 
contamination, which arc estimated to be < 10%. 

The redshift selection function of the X-ray sources, 
derived by using the soft-band luminosity function of 
Ebrero et al. (2009b) which takes into account the re- 
alistic luminosity dependent density evolution of the X- 
rays sources, predicts a characteristic depth of z ~ 1. 

In Figure 1, we present the X-ray AGN angular cor- 
relation function of the Ebrero et al. (2009a) analysis. 
The solid points corresponds to the observed angular 
correlation function while the solid line represents the 
theoretical angular correlation function for the best fit- 
ting cosmological model (see further below) . The insert 
panel of Fig.l shows the residual, Aw(^), between obser- 
vations and theory. There appears to be an interesting 
sinusoidal variation with 9, which merits further inves- 
tigation. Unaccounted non-linear effects, at the smallest 
angular separations, could be the cause of the large Aw 
values at ^ < 80 . 



2.2 From Angular to Spatial Clustering 

We briefly present the main points of the method used 
to put cosmological constraints using the angular clus- 
tering of some extragalactic mass-tracer. A first impor- 
tant important ingredient is the use of Limber's formula 
which relates the angular, w{9), and the spatial, £,(r), 
correlation functions. Assuming flatness the Limber's 
equation can be written as: 



w{e) = 2 
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(1) 



where 4>{x) is the AGN redshift selection function (the 
probability that a source at a distance x is detected 
in the survey) and x is the coordinate distance re- 
lated to the redshift through x{z) = c j^^ '^{y) ^it^ 
H[z) = i?„[n„,(l + zf + Q.q{1 + ^)3{i+w)]i/2 
Q,Q = 1 — Q,m- Also, r is the physical separation be- 
tween two sources having an angular separation, 9, 
which in the small angle approximation is given by 
r ~ (1 z)-^ [u^ + x^e'^y^ (with u the line-of-sight 
separation of any two sources). The number of AGN 
within a shell {z, z -\- dz) is given by: 
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where (5a;s(~ 125.5 deg^ is the effective solid angle of 
the survey, E{z) = H{z)/Ho and is the comoving 
AGN number density at « = 0. The source redshift 
distribution AN/Az, as already mentioned previously, is 
estimated by integrating the appropriate Ebrero et al. 
(2009b) luminosity function, folding in the area curve 
of the survey. 

Inserting eq.(2) into eq.(l), we have after some al- 
gebra that: 



w(9) 



^ c \N Az 



05 




e(arcsec) 

Figure 1. The two point angular correlation function of the 
soft band (Ebrero et al. 2009a). The solid line reijrcsents 
the theoretical angular correlation function for the best fit- 
ting {Clm = 0.26 and w = —0.93) cosmological model. Insert 
Panel: The residual angular correlation function between ob- 
servations and theory with 2cr uncertainties. 



The spatial correlation function can be written as: 
C{r,z) = (1 + z)-^^+'^b''{z)CoM{r), with CoM^r) in- 
dicating the predicted spatial correlation function of 
the underlying matter distribution (see below) and e 
parametrizing the type of AGN clustering evolution 
(eg. de Zotti et al. 1990) and following Kiindic (1997) 
and Basilakos & Plionis (2005; 2006) we use here the 
constant in comoving coordinates clustering model, ie., 
€ = —1.2. Also, b{z) is the evolution of the linear bias 
factor which is an essential ingredient for cold dark 
matter (CDM) models in order to reproduce the ob- 
served mass-tracer distribution (cf. Kaiser 1984; Davis 
et al. 1985; Bardeen et al. 1986; Benson et al. 2000). 
In the current analysis we use our bias evolution model 
(Basilakos, Plionis & Ragone-Figueroa 2008), which is 
based on the solution of a second order differential equa- 
tion derived by using linear perturbation theory and the 
Priedmann-Lemaitre solutions of the cosmological field 
equations. Our model was initially presented in Basi- 
lakos & Plionis (2001; 2003) and has been recently ex- 
tended to include the effects of halo interactions and 
merging (for details see Basilakos et al. 2008). 



2.3 Theoretically Predicted Clustering 

We estimate the theoretically predicted spatial corre- 
lation function of the underlying matter distribution, 
^DM('r), from the Fourier transform of the spatial power 
spectrum P{k): 



$dm( 
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(4) 



E{z)Az / ^{r,z)Au . (3) 



where P{k) denotes the power of the matter fluctuations 
linearly extrapolated to the present epoch. We consider 
the CDM power spectrum, P{k) = Pofe"T^(fc), with 
T{k) the CDM transfer function and and n ~ 0.96 fol- 
lowing the 5-year WMAP results (Komatsu et al. 2009). 
In order to deflne the functional form of the power spec- 
trum, we utilize the transfer function parameterization 
as in Bardeen et al. (1986), with the approximate cor- 
rections given by Sugiyama (1995). The rms fluctuations 
of the linear density field on mass scale M is: 
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a{M) = 
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(5) 



where 

the window function is given by: W{kR) = 3(sinfc_R — 
kRcoskR) / [kRf and R = (3Af/47rpo)^^'^ The parame- 
ter po denotes the mean matter density of the universe 
at the present time (po = 2.78 x lQ^^Q.mh? M^Mpc'^). 
The normahzation of the power spectrum is given by: 
Po = 27rV| [/J"r2(fc)fc"+2Vl/2(fci?8)dfc]"^ where as is 
the rms mass fluctuation on Rs = 8h~^ Mpc scales and 
for which we use the WMAP5 value of ere — 0.8 (Ko- 
matsu et al. 2009). It is worth noting that we also use the 
non-linear corrections introduced by Peacock & Dodds 
(1994). 



2.4 Cosmological Constraints 

In order to constrain the cosmological parameters we 
use, as in Basilakos & Plionis (2005), a standard like- 
lihood procedure and compare the measured XMM soft 
source angular correlation function (Ebrero et al. 2009a) 
with the predictions of different spatially flat cosmologi- 
cal models. To this end we use the likelihood estimator*, 
defined as: £agn(p) oc exp[-xiGN(p)/2] with: 



2 

Xagn 



(p) = E 



[Wth{Oi,p) - Wobs{9i)f 



(6) 



where p is a vector containing the cosmological param- 
eters that we want to estimate, ai is the uncertainty of 
the observed angular correlation function and ae, cor- 
responds to the width of the angular separation bins. 

As we have previously mentioned, we work within 
the framework of a flat cosmology with primordial adi- 
abatic fluctuations and baryonic density of Qhh^ = 
0.022(±0.002) (eg. Komatsu et al. 2009), while utilizing 
the HST key project results of Freedman et al. (2001) we 
flx the Hubble constant to ~ 0.71 (see also Komatsu 
et al. 2009). Note that since we flx, in the following anal- 
ysis, the values of both h and Qh, we do not take into 
account their quite small uncertainties. 

The corresponding statistical vector that we have 
to flt is: p = (ftm,w, Mh), where Mh is the host dark 
matter halo mass, which enters in our biasing evolu- 
tion scheme (see Basilakos, Plionis & Ragone-Figueroa 
2008). Note also that the normalization of the power 
spectrum, as, could have been left as a free param- 
eter (see Basilakos & Plionis 2006), but in this work 
we choose to use the well established WMAP5 value 
(see previous section). We sample the various param- 
eters as follows: the matter density ft,n G [0.01, 1] in 
steps of 0.01; the equation of state parameter w G 
[—1.6, —0.34] in steps of 0.01 and the parent dark matter 
halo Mfe/lO^^/r^Mo £ [0.1, 4] in steps of 0.01. 

We find that the likelihood function of the soft X- 
ray sources peaks at f2„i — 0.26 ± 0.05, w = —0.931" jg 
and Mh = 2to't x 10^^ H'^Mq, with a reduced of 



* Likelihoods are normalized to their maximum values. 



~ 4. Such a large value is caused by the small ^(6^) un- 
certainties in combination with the observed modula- 
tion (see insert panel of Fig.l). Had we used a 2a 'w{6) 
uncertainty in eq.(6) we would have obtained roughly 
the same constraints and a reduced of ~ 1 (see up- 
per right panel of Fig.2). 

In the upper-left panel of Figure 2 we present the 
current constraints in the (fim, w) plane by marginaliz- 
ing our solution over Mh (thick solid lines). For com- 
parison reasons we also show our previous solutions of 
Basilakos & Plionis (2005; 2006), which where based on 
the shallow XMM/2dF survey (~ 2.3 deg^) which con- 
tains only 432 point sources (with an effective flux-limit 
of fx > 2.7 X 10"^" erg cm"^ s"^). In particular, the 
dotted lines correspond to a solution using erg — 0.93 
(Basilakos & Plionis 2005), while the dashed-lines to 
the corresponding solution for as ^ 0.74 (Basilakos & 
Plionis 2006). Comparing our present analysis with our 
previous results it becomes evident that with the cur- 
rent high-precision X-ray AGN correlation function of 
Ebrero et al. (2009a) we have achieved to place simulta- 
neously quite stringent constraints on both, w and Qm- 

Regarding other analyses of cosmological data, it 
is interesting to note that Davis et al. (2007) using 
the combined analysis of the SNIa-|-BAO-|-CMB found 
Qm = 0.27 ± 0.04 and w = -1.01 ± 0.15, while a sim- 
ilar analysis of Kowalski et al. (2008), using a newer 
SNIa compilation, provided flm = 0.274 ± 0.016 with 
w= —0.969 ±0.06 (see also corresponding results in Ko- 
matsu et al. 2009 and Hicken et al. 2009). We would 
like to stress that despite the fact that we use a single 
cosmologically relevant experiment, ie., the observed an- 
gular correlation function of the soft X-ray sources, our 
results coincide within la with the results of the joint- 
analysis, discussed previously. Therefore, the X-ray se- 
lected AGNs appear to be ideal tools for extracting cos- 
mological information. In order to further illustrate such 
a claim, we perform further below a direct comparison 
between our results with those derived by other single 
cosmological datar-sets. 

2.5 Comparison with other Cosmological Data 

We present here a comparison between the {Qrn,'w) 
solution-space provided by our analysis of the XMM- 
Newton X-ray sources with those derived using other 
cosmological data. The goal is to give the reader the 
opportunity to appreciate the relative strength and pre- 
cision of the different methods. Therefore, we present 
in Figure 2 the la, 2a and 3a confldence levels in the 
(fimjw) plane, provided by the different cosmological 
data presently available. These are: (a) the Hubble re- 
lation based on the latest sample of 397 supernovae of 
Hicken et al. (2009); (b) the dimensionless distance to 
the surface of the last scattering R = 1.71 ± 0.019 (Ko- 
matsu et al. 2009); (e) the baryon acoustic oscillation 
distance aX, z = 0.35, A = 0.469±0.017 (Eisenstein et al. 
2005; Padmanabhan, et al. 2007); (d) the Hubble func- 
tion (Hz) derived directly from 9 early type galaxies at 
high redshifts (Simon, Verde & Jimenez 2005) and (e) 
the Hubble relation based on a sample of 69 Gammarray 
bursts (Cardone, Capozziello & Dainotti 2009). It is ev- 
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Figure 2. Likelihood contours (Icr, 2(t and 3o") in the 
(n,„ , w) plane. The upper left panel shows the results based 
on the current (solid lines) and previous XMM source clus- 
tering analyses. In the upper right panel we show the results 
based on the SNIa (solid lines) and the XMM clustering re- 
sults and a 2(7 clustering uncertainty (dashed lines) . The rest 
of the panels correspond to the cosmological data indicated. 
The contours are plotted where — 21ni2/i2max is equal to 2.30, 
6.16 and 11.83. 



idcnt that the X-ray AGN clustering likchhood analysis 
alone puts the most stringent constraints on the value 
of the equation of state parameter. To be fair however 
we must stress that in our analysis we have a priori im- 
posed the value of ag to that of the WMAP5 (Komatsu 
et al. 2009), which docs not enter the other cosmologi- 
cal tests. However, the size of the solution space, for any 
plausible value of as, is as small as for the nominal case 
used, the only difference would be a tilt of the contours, 
as can be appreciated by the two XMM/2dF solution 
spaces in the upper left panel of Fig. 2. 

Note, that if we increase the uncertainty of the ob- 
served X-ray source angular correlation function by a 
factor of 2, then the resulting contours (thick dashed 
lines in the right panel of figure 2) match closely those 
of the most recent SNIa analysis. In a forthcoming paper 
(Plionis et al. in preparation) we will present details of 
a joint-analysis of our X-ray selected AGN results with 
that of all other cosmologically relevant data. 

Lastly, we have to caution the reader of two (rea- 
sonable, we believe) assumptions that enter a priori in 
our analysis: (1) that the clustering evolution of X-ray 
selected AGN is constant in comoving coordinates (the 
effects of other evolution models has been investigated 
in Basilakos & Plionis 2005) and (2) that the Basilakos 
et al. (2008) bias evolution model is the appropriate one, 
which is supported by a comparison with N-body sim- 



ulations and available clustering data (see figures 1 and 
3 of Basilakos et al. 2008). 



3 CONCLUSIONS 

We have utilized the recent determination of the cluster- 
ing properties of high-z X-ray selected AGN, identified 
as soft (0.5-2 keV) point sources (Ebrero et al. 2009a), in 
order to constraint the main cosmological parameters. 
We find that the X-ray AGN clustering likelihood anal- 
ysis alone, within the context of flat cosmological mod- 
els, can place tight constrains on the main cosmological 
parameters (f2m,w), and indeed relatively tighter than 
any other single observational method todatc. The cur- 
rent analysis provides a best spatially fiat model with 
Qm = 0.26 ± 0.05 and w= -0.93lo;l9. 
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